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We determine abundances from the absorption spectrum of the magnetic Herbig Ae star HD 190073 (V1295 Aql). The 
observations are primarily from HARPS spectra obtained at a single epoch. We accept arguments that the presence of 
numerous emission lines does not vitiate a classical abundance analysis, though it likely reduces the achievable accuracy. 
Most abundances are closely solar, but several elements show departures of a factor of two to three, as an earlier study has 
also shown. We present quantitative measurements of more than 60 emission lines, peak intensities, equivalent widths, and 
FWHM's. The latter range from over 200 kms _1 (Ha, He D3) down to 10-20 km s _1 (forbidden lines). Metallic emission 
lines have intermediate widths. We eschew modeling, and content ourselves with a presentation of the observations a 
successful model must explain. Low-excitation features such as the Na I D-lines and [O I] appear with He I D3, suggesting 
proximate regions with widely differing T c and N c as found in the solar chromosphere. The [O I] and [Ca II] lines show 
sharp, violet-shifted features. Additionally, [Fe II] lines appear to be weakly present in emission. 
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1 Introduction 

The intriguing spectrum of the magnetic Herbig Ae star 
HD 190073 (V1295 Aql) has attracted the attention of clas- 
sical as well as modern spectroscopists (Merrill 1933; 
Swings & Struve 1940; Catala et al. 2007, henceforth CAT). 
Pogodin, Franco & Lopes (2005, henceforth P05) give a de- 
tailed description of the spectrum along with a historical re- 
sume of investigations from the 1930's. We note that the 
nature of HD 190073 as a young, Herbig Ae star became 
widely recognized some three decades after Herbig's (1960) 
seminal paper. 

HD 190073 was included among the 24 young stars 
studied for abundances by Acke & Waelkens (2004, hence- 
forth, AW). In this important paper, the authors made the 
bold assumption that abundances could be determined for 
stars of this nature using standard techniques-plane paral- 
lel, one dimensional models, in hydrostatic equilibrium. The 
models were used to obtain abundances from absorption 
lines with equivalent widths less than 150 mA. These as- 
sumptions might very well be questioned. Material is being 
accreted by these young stars, and the infall velocities are 
thought to be near free-fall, several hundred kms -1 . Does 
this infall produce shocks and heating of the atmospheres 
that could invalidate models that neglect such complica- 
tions? 

AW nevertheless proceeded. Although they did not state 
this explicitly, the justification for their assumptions is em- 
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pirical, and may be found in their results. Basically, these 
are the fact that their approach yields entirely reasonable 
stellar parameters and abundances, including agreement 
from lines of neutral and ionized elements. Stated simply, 
their assumptions led them to self-consistent results. We 
make these same assumptions in the present work, taking 
some comfort in the fact that self consistency is all one ever 
has in science. 

While AW's studies were both competent and thorough, 
better observational material is currently available, making 
it possible to use systematically weaker lines, and to study 
more elements. We have also made use of the wings of the 
Balmer lines, not used by AW. 

The lower Balmer lines have central emissions. In the 
case of Ha, the emission dominates the feature. The Balmer 
lines and especially Ha, have been extensively studied (e.g. 
P05; Cuttela & Ringuelet 1990). 

In the present paper we also study the weaker metallic 
emission lines, to provide information on the physical con- 
ditions where this emission occurs. This was discussed by 
CAT who were primarily concerned with the magnetic field 
of HD 190073. They also give a detailed qualitative descrip- 
tion of the metallic emission lines (primarily Ti II, Fe i/ll (cf. 
Sect. [6]and following). 

Hubrig et al. (2006, 2009) reported a longitudinal mag- 
netic field of 84±30 Gauss, up to 104±19 Gauss, while 
CAT found a longitudinal field of 74±10 Gauss. 
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Fig.l The HARPS spectrum (ADR HARPS. 2008-1 1- 
10T23:43:14.386^_SID_A) and averaged spectra in the re- 
gion of the Mg II doublet A4481 A. The HARPS (gray and 
red in online version) spectrum has been displaced slightly 
downward for display purposes. 
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Fig. 2 The region from AA4000 to 4050 A shows numer- 
ous measurable absorption features. Note the broad absorp- 
tion at A4026 A, which is He I. Fe I A4045 A shows strong 
emission as well as absorption. 

2 Observations 

We downloaded 8 HARPS spectra from the ESO archive, 
all obtained on 1 1 November 2008 within 74 minutes of 
one another. These were averaged, binned to 0.02 A and 
mildly Fourier filtered. The resulting spectrum has a signal- 
to-noise ratio of 350 or more. The resolution of HARPS 
spectra usually cited as over 100000, is not significantly 
modified for our purposes on the averaged spectrum, as 
Fig- [U illustrates. 

UVES spectra, obtained on 18 September 2005 cover 
the region from 3044 to 10257 A. They were used for spe- 
cial features (e.g. the [Ca II] lines), but not for abundances. 

3 Reduction 

The averaged HARPS spectrum was measured for 1796 
lines. The UVES spectrum was also measured for line iden- 
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tifications in the region AA3054-3867 A. We measured 760 
absorption lines, which were often severely affected by 
emission. 

Many absorption lines, especially weak ones, were not 
significantly affected by emission, and suitable for abun- 
dance determinations. Figure |2] shows a typical region with 
many relatively unperturbed absorptions. Preliminary, au- 
tomated identifications were made, and wavelength coin- 
cidence statistics (WCS, Cowley & Hensberge 1981) were 
performed. A few spectra not investigated by AW were an- 
alyzed: He I, Na I, Al I, Si I, S I, S II, Co I, Mn I, Mn II, 
Ni I, Zn I, and Sr II. We found no exotic elements, such as 
lanthanides, or unusual 4d or 5d elements. 

Lines were chosen for equivalent width measurement 
with the help of the automated identification list, which lists 
plausible identifications within 0.03 A of any measured fea- 
ture. Blends were rejected. Usually, we avoided lines with 
equivalent widths greater than 20 mA but in order to com- 
pare our measurements with those of AW, we included a few 
stronger lines. 

A comparison of measurements is given in Fig. [3] Gen- 
erally, the measurements agree well with one another, and 
differences can usually be explained by judgments of where 
to draw the continuum when a line is partially in emission, 
or there is emission close by. Differences in the case of 
one of the solid circles is surely due to emission, as Ti II 
A4398 A falls between two strong emission lines. The other 
solid point is for O I A3947 A. This is apparently a misiden- 
tification. Note that Fig.[3]is not logarithmic. 

4 The model atmosphere and abundance 
methods 

The methods used to obtain abundances from the equivalent 
widths, including model atmosphere construction are ex- 
plained in some detail in two previous papers (Cowley et al. 
2010a,b). Briefly, the T{t 5000 ) from Atlas 9 (Kurucz 1993) 
as implemented by Sbordone et al. (2004) was used with 
Michigan software to product depth-dependent models. The 



© 20 1 1 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (201 1) 



3 



Table 1 Abundances in HD 190073 from the current 
study and AW. 



Ion 
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3 


-9.86 


-9.72 



effective temperature and gravity were selected from ioniza- 
tion and excitation equilibrium as well as fits to the wings 
of U/3-US. 

We have adopted a somewhat lower temperature than 
used by AW, 8750 K, and log 5 = 3.0. The former used 
9250 K, and logy = 3.5. We also used a lower microturbu- 
lence, 2 km s" 1 , compared to AW's 3 km s -1 , but this is not 
important for most of our weaker lines. Oscillator strengths 
were taken from the modern literature when possible, or 
from compilations by NIST (Ralchenko 2010, preferred) or 
VALD (Kupka et al. 1999). Default damping constants were 
used as in the studies cited, but they are unimportant for 
weak lines. 

5 Abundances 

The AW abundances have been converted from differential 
values, using Anders & Grevesse (1989) abundances, which 
AG adopted. Our abundances (see Table [TJ refer to the As- 
plund et al. (2009) scale. The case is not strong that any of 
these abundances differ significantly from the solar abun- 
dance. Nevertheless, we have highlighted in bold face some 
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Fig. 4 He I D 3 line in the spectrum of HD 190073. The 
central absorption is arguably photospheric, and agrees in 
shape and strength with a calculated absorption profile. 

elements that deserve additional attention. Nitrogen, in par- 
ticular, deserves attention, as it has been found in excess in 
the Herbig Ae star HD 101412 (Cowley et al. 2010a). As- 
terisks mark cases where AW and the present work agree 
on possibly significant departures from solar abundances. 
NLTE effects could also be responsible for some non-solar 
abundances (Kamp et al. 2001). 

5.1 Neutral helium 

The helium abundance is from AA4026 (see Fig. [2]) and 
4713 A. Both lines are weak, but in excellent agreement 
with one another. However A447 1 A was also found in ab- 
sorption, and analyzed. The value of log(He/iVtot) from 
this line was found to be —1.53, some 0.4 dex below the 
mean of the other two lines. We have chosen to disregard 
this value, as possibly weakened by partial emission. Should 
it be included in the average, we find log(He/iV tot ) = 
—1.29 ± 0.21 (sd), still solar, within the errors. The D 3 line 
(A5876 A) of He I is in emission, and included in Table [2] 
It was observed at numerous phases by P05, whose obser- 
vations (see their Fig. 3) do not show the central rever- 
sal clearly seen in our Fig. [4] This feature was measured 
at A*5875.60, virtually unshifted from the expected photo- 
spheric position A5875.64. Moreover, an LTE synthesis us- 
ing Voigt profiles and an assumed solar abundance matches 
the observed absorption in shape and strength. Given the 
likelihood of NLTE, it is unclear how seriously to take this 
result. Nevertheless, the D3 absorption is consistent with 
photospheric absorption, and a solar abundance. 

The D3 emission is remarkably similar in morphology to 
the metallic emissions (with unshifted absorptions), though 
it is much broader. It resembles the P05 illustrations at 
phases 'a' or 'e', with a maximum shifted somewhat to the 
red, and a longer violet tail. 
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Fig. 5 The emission/absorption, metallic-line spectrum of 
two Herbig Ae stars, contrasted. Note the proclivity of the 
intrinsically stronger lines to be in emission in HD 190073. 

6 The emission spectrum 

A second focus of the current paper is the emission spec- 
trum, in particular, permitted and forbidden metallic lines. 

Previous studies (P05, CAT) provide detailed descrip- 
tions of numerous metallic emission lines, primarily of Fe I 
and II, Ti II, and Cr II, as well as the Na D-lines, and have 
illustrations similar to the upper spectrum of Fig. [5] As the 
P05 work has a temporal dimension lacking in the present 
study, we briefly summarize their findings. The emissions 
show mild temporal variations both in strengths and widths. 
The profiles are somewhat asymmetric, with their peak in- 
tensities generally very slightly red shifted with respect to 
the photospheric absorption spectrum. The widths of the 
features are significantly larger than these shifts, and a con- 
siderable fraction of the emission is shifted to the blue. 

The P05 observations are all in good agreement with 
the current findings, which we take to be a representative 
sample. Table [2] gives measurements of the peak intensi- 
ties, equivalent widths, and FWHM for the strong relatively 
unblended emission lines as they appear on our averaged 
HARPS spectra. Multiplet numbers follow the spectrum 
designation. The intensities are in units of the continuum, 
and the equivalent widths are the areas of the emissions 
above the continuum, which is assumed to have unit inten- 
sity. The measurements are from segments fitted by eye to 
the emission lines, as illustrated in Fig. [6] 

All lines in the table had central absorptions. In measur- 
ing the /"-values, an attempt was made to interpolate over 
this absorption, so the 7°-value is a few per cent higher than 
the maximum of the emission. The accuracy of the mea- 
surements vary. Repeated measurements show consistency 
for FWHM are generally within 10 %. Underlying emission 
from other lines is the prime cause of the uncertainty. 

6.1 Resume: the permitted emissions 

We summarize salient properties of the permitted emission 
lines: 
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Table 2 Maximum intensity measurements, equivalent 
widths, and full widths at half maximum for selected emis- 
sion lines. 
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Fig. 6 Segment fits (black) to the emission from Mg 1-2 
A5183.60 A, the strongest of the Mg lb lines. Observations: 
gray (red online) with dots. The portion of the fitted curve 
near the central absorption is a by-eye estimate of the miss- 
ing part of the profile. The maximum designated 1°, and 
other properties of similar interpolations are given in Ta- 
ble0 

- The centers of gravity are shifted by ca. 5 km s -1 to the 
red. 

- The profiles are somewhat asymmetric, with a longer 
violet than red tail. 

- The central absorption wavelengths are photospheric 
within the errors of measurement. That is, the radial 
velocities of the weaker photospheric absorptions agree 
with those of the central absorption components of the 
emission lines. 

- The emission lines are the intrinsically strongest lines. 
Weaker lines show weaker emission, until, for the weak- 
est lines, the observed features are all in absorption. 

CAT suggested the emissions arose in conditions similar to 
those of a photosphere. They suggested a heated region with 
densities and temperatures in the range of 10 13 -10 14 cm -3 
and 15 000-20 000K. 

Their value of 65 kms -1 as a typical FWHM for the 
emissions agrees well with our measurements (Table [2]). 
The origin of this velocity, however, is not readily apparent. 
They speculate that these velocities are due to a supersonic 
turbulence. Such "turbulence" might arise from the roil of 
accreting material settling onto the photosphere. 

6.2 Forbidden lines 
6.2.1 The [O i] Unes 

Both AA6300 and 6363 [O I] are present as well as the au- 
roral transition A5577 A. In addition to what we shall call 
the main features, all three lines show faint, sharp, "satel- 
lite" components shifted to the violet by ca. 25 kms -1 . This 
structure is illustrated in Figs. [7] and |T0j The main [O I] 
features are roughly one third the width of the typical per- 
mitted metallic-line features, but their peak intensities have 
comparably small red shifts of ca. 5 kms -1 . It is plausible 



x 

1) 
> 
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Wavelength [A] 

Fig. 7 The nebular ( 3 P 2 - 1 D 2 ) [O I] transition A6300 
A. The long arrow points to the laboratory position at 
6300.30 A. The stellar feature seems slightly red shifted. 
The short arrow points to a narrow, blue-shifted satellite 
feature that is also seen at the same displacement in [O I] 
AA6363 and 5577 A. The strong absorption lines are atmo- 
spheric. 
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Fig. 8 The [Ca II] lines AA729 1 (black) and 7324 A (gray, 
red in online version) from multiplet 1-F. Wavelengths for 
A7324 A are at the top abscissa. The arrow indicates a pos- 
sible narrow component of A7291 A corresponding to those 
seen in [O I]. 

to assume the [O I] arises in a region further from the star, 
and therefore of lower density than the gas giving rise to the 
permitted metallic lines. 

The satellite emissions may arise in a polar stream, if we 
assume the system is viewed pole on. The velocity, however, 
is not high. The satellite features of all three [O I] lines is 
ca. —25 kms -1 . 

6.2.2 The [Ca n-Fl] doublet 

Hamann (1994) has noted the presence of [Ca II] in a num- 
ber of young stars, including the Herbig Ae V380Ori. We 
are not aware that the lines have been previously noted in 
HD 190073. They are also seen in supernovae (Kirshner & 
Kwan 1975) and extragalactic spectra. (Donahue & Voit 
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Fig. 9 Forbidden Fe II lines from multiplet 19-F. The 
wavelength scale for the weaker of the two lines, A5376 A 
(gray, red in online version) is at the top of the figure. Ver- 
tical lines mark the positions of wavelengths derived from 
the atomic energy levels. 

1993). Merrill (1943) reported [Ca II] lines in emission in 
the peculiar hydrogen-poor binary u,Sgr (see also Green- 
stein & Merrill 1946). 

We find definite emissions at the positions of the for- 
bidden 2 Si/2 - 2 D3/2,5/2 transitions. The air wavelengths, 
determined from the energy levels, are 7291.47 and 
7323.89 A. These features were identified on a the UVES 
spectrogram. The spectrum (Fig. [8J is too noisy or blended 
to show the presence of satellite features of A7324 A, but it 
might be present for A7291 A. Measurements of the main 
features are included in Table [2] The maxima are shifted 
to the red by 2 to 4 kms^ 1 , in general agreement with the 
[Ol] and metallic lines. The FWHM agree with those of 
other forbidden lines. 

6.2.3 [Fe II] 

Several workers have discussed [Fe II] emission lines in 
Herbig Ae/Be stars (Finkenzeller 1985; Donati et al. 1997). 
We find a definite, sharp emission feature with a maximum 
measured at A*5 158.84 A. This wavelength is close to that 
of [Fe II-19F], A5 158.78 A. (Laboratory positions for [Fe II] 
are from Fuhr & Wiese (2006) rather than the RMT). An- 
other line from this multiplet is weakly present (Fig. [9]). 
Both features are seen on the unaveraged HARPS spectra 
as well as UVES spectra taken some 3 years previously (see 
Sect. [2}. Several other lines in Multiplet 19-F are arguably 
present (AA5261, 5296, 5072 A), other lines are masked by 
blends or in a HARPS order gap. The line A5158 A is en- 
tered in Table [2] It has a FWHM comparable to that of the 
other forbidden lines. 

We found no other [Fe II] features that could be said 
to be unambiguously present. The well-known [Fe II] line 
A4244 A is at best, marginally present. 
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Fig. 10 The auroral ( 1 D 2 - 1 S ) [O I] transition A5577 
A (gray, red in online version). The long arrow points to 
5577.34 A, the NIST wavelength. The maximum, and center 
of gravity of the line shifted slightly to the red, as are the 
nebular lines (Pig. [7J. The shorter arrow points to the sharp, 
blue-shifted feature with the same displacement as the sharp 
component of [O i] AA6300 and 6363 A. 

6.2.4 Physical conditions from the forbidden lines 

"Critical electron densities" are obtained from observed for- 
bidden transitions by equating the Einstein spontaneous de- 
cay rate to the collisional deexcitation rate. Typical values 
are given by Draine (201 1) for T = 10 000K. For the [O I] 
^a-level (the upper level of AA6300 and 6363 A), the crit- 
ical N e is 1.6xl0 6 cm~ 3 . We calculate a similar critical 
density for the [Ca II] lines with the help of rates calculated 
by Burgess et al. (1995), and a lifetime of the 2 D term given 
by NIST (Ralchenko et al. 2010). 

We see no evidence of the [S II-2F] pair at AA6717 and 
673 1 A, though they have been observed in Herbig Ae/Be 
stars (Corcoran & Ray 1997). For this pair, Draine gives 
critical densities of 10 3 -10 4 cm -3 . We conclude the for- 
bidden lines we do see arise in a region where the electron 
density is between 10 4 and 10 6 ~ 7 cm~ 3 . 

When the nebular as well as auroral transitions of [O I] 
are available, the ratio may allow one to determine values of 
T e and 7V C compatible with the observation. The average of 
three measurements on A5577 gives W = 0.0042 A, which 
with W = 0.056 A for A6300 yields a ratio of 13.3. If we 
assume the excited levels of O Iarise from electron excita- 
tion, we may interpolate in the plot of Gorti et al. (201 1) to 
find acceptable the values given in Table[3] With the electron 
density constraint given above, we find temperatures in the 
range 7500 to 10000 K for the volume where the forbidden 
lines are formed. 

An alternate interpretation of [O I] emission in Her- 
big Ae/Be systems is discussed by Acke, van den An- 
cker & Dullemond (2005). In their model, the excited O 
ilevels arise primarily from the photodissociation of the OH 
molecule. 
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Table 3 Values of T c (in K) and N c (in cm 3 ) compatible 
with the observed A6300/5577-ratio = 13.3. 





logiV c 


5000 


8 


7500 


7 


10000 


6.5 


12000 


6.2 



Acknowledgements. We are grateful for the availability of the ESO 
archive. This research has made use of the SIMBAD data base, 
operated at CDS, Strasbourg, France. Our calculations have made 
extensive use of the VALD atomic data base (Kupka et al. 1999), 
as well as the facilities provided by NIST (Ralchenko et al. 2010). 
CRC thanks colleagues at Michigan for many helpful suggestions. 
Jesus Hernandez suggested that we examine the forbidden lines. 

References 

Acke, B., Waelkens, C: 2004, A&A 427, 1017 (AW) 
Acke, B., van den Ancker, M.E., Dullemond, CP.: 2005, A&A 
436, 209 

Anders, E., Grevesse, N.: 1989, Geochim. Cosmochim. Acta 53, 
197 

Asplund, M., Grevesse, N., Sauval, A.J., Scott, P.: 2009, ARA&A 
47, 481 

Burgess, A., Chidichimo, M.C., Tully, J.A.: 1995, A&A 300, 627 
Catala, C, Alecian, E., Donati, J.-E, et al.: 2007, A&A 462, 293 
(CAT) 

Corcoran, M., Ray, TP: 1997, A&A 321, 189 
Cowley, C.R., Hensberge, H.: 1981, ApJ 244, 252 
Cowley, C.R., Hubrig, S., Gonzalez, J.F., Savanov, I.: 2010a, A&A 
523, 65 



Cowley, C.R., Hubrig, S., Palmeri, P., et al.: 2010b, MNRAS 405, 
1271 

Cuttela, M., Ringuelet, A.E.: 1990, MNRAS 246, 20 

Donahue, M., Voit, G.M.: 1993, ApJ 414, L17 

Donati, J.-E, Semel, M., Carter, B.D., Rees, D.E., Cameron, A.C.: 

1997, MNRAS 291, 658 
Drain, B.T.: 2011, Physics of the Interstellar and Intergalactic 

Medium, Princeton Univ. Press, Princeton, NJ 
Finkenzeller, U.: 1985, A&A 151, 340 

Fuhr, J. R. & Wiese, W. L.: 2006, J. Phys. Chem. Ref. Data, 35, 
1669 

Greenstein, J.L., Merrill, P.W.: 1946, ApJ 104, 177 
Hamann, E: 1994, ApJS 93, 485 
Herbig, G.H.: 1960, ApJS 4, 337 

Hubrig, S., Yudin, R.V., Scholler, M., Pogodin, M.A.: 2006, A&A 
446, 1089 

Hubrig, S., Stelzer, B„ Scholler, M., et al: 2009, A&A 502, 283 
Kamp, I., Iliev, I.Kh., Paunzen, E., et al.: 2001, A&A 375, 899 
Kirshner, R.P, Kwan, J.: 1975, ApJ 197, 415 
Kupka, E, Piskunov, N.E., Ryabchikova, T.A., Stempels, H.C., 

Weiss, W.W.: 1999, A&AS 138, 119 
Kurucz, R.L.: 1993, ATLAS9 Stellar Atmosphere Programs and 

2 km/s grid, CD-Rom, No. 13, Smithsonian Ap. Obs., Cam- 

bridge/MA 
Merrill, P.W.: 1933, ApJ 77, 51 
Merrill, P.: 1943, PASP 55, 242 

Pogodin, M.A., Franco, G.A.P, Lopes, D.E: 2005, A&A 438, 239 

Ralchenko, Yu., Kramida, A.E., Reader, J., NIST ASD team: 
2010, NIST Atomic Spectra Database, version 3.1.5, 
http://physics.nist.gov/as of [2010, Feb 10], National Institute 
of Standards and Technology, Gaithersburg, MD 

Sbordone, L., Bonifacio, P., Castelli, E, Kurucz, R.L.: 2004, Mem. 
Soc. Astron. Ital. Suppl. 5, 93 

Swings, P., Struve, O.: 1940, ApJ 91, 546 



www.an-journal.org 



© 201 1 WILEY-VCH Verkg GmbH & Co. KGaA, Weinheim 



